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Propagation of acoustic waves within a continuum of inviscid, compressible, heat-conducting 
fluid is evaluated in detail in terms of both frequency (steady-state) and time dependent (transient) 
functions. The analysis reveals that when the value of the ratio of specific heats, y, lies between one 
and two, the apparent steady-state solutions are conjugate to unstable, or regenerative, "transient" 
solutions, and, thus, are unacceptable. Propagation is stable for other values (y = l,y 3= 2). The common 
assumption that the steady-state phase velocity varies continuously with increasing frequency from 
adiabatic to isothermal values is shown to be invalid, except when y=2. 
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1 . Introduction 

This paper presents an analysis of the propagation of acoustic waves of small amplitude within 
a continuum of inviscid, compressible, heat-conducting fluid. The analysis gives the conditions 
imposed on such a fluid for the propagation to be stable. In this usage "stable" propagation is that 
manifested by a "stable" fluid, which has no internal energy sources (i.e., a passive fluid). With 
acoustic propagation the concept of stable propagation implies that for a medium subjected to 
sinusoidally time-dependent excitation commencing at some time, the response at any point in 
the medium will approach a bounded steady state. Although the requirement for stable propagation 
is obvious for any physically realistic fluid (as far as we know), its full significance has apparently 
been overlooked. In general the inclusion of this requirement will impose a relationship [1J 1 between 
certain thermophysical properties of the fluid. 

In an earlier paper [1] stability criteria, given as inequalities in terms of certain thermophysical 
properties, were developed for viscous, heat-conducting fluids, which are a more general class. 
These criteria were developed without directly analyzing the explicit time dependence of the exact 
solution. In the following presentation the inviscid fluid, which is simpler to analyze, is examined 
in detail with respect to both steady-state and transient functions. The analysis reveals some 
contradictions in present generally accepted concepts with respect to propagation in inviscid 
fluids. Propagation in these fluids has been assumed to be stable and the steady-state propagation 
functions have been assumed to vary continuously from adiabatic conditions at low frequencies to 
isothermal at high frequencies. Our analysis reveals that propagation is unstable when the value 
of the ratio of specific heats, y, lies between one and two, and continuous variation in the steady- 
state propagation functions from adiabatic conditions at low frequencies to isothermal at high 
frequencies is not possible in any stable fluid, except when y = 2. The unstable class of fluids may, 
of course, be stabilized [1] by the inclusion of viscosity of sufficient magnitude within the model. 

Considerable discussion on propagation in heat-conducting fluids has appeared in the literature 
over the past 100 years, for which some of the more recent examples are given in references [2— 9J. 



1 Figures in brackets indicate the literature references at the end of this paper. 
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List of Symbols 

a = pcgCv/K = C\IkB s 
&=pcJ/(t/+4r>/3) 

B* = B' — iB" = complex compressibility 
B' = Re B* 
B"=-ImS* 

5.s = adiabatic (isentropic) compressibility 
Br — yB s — isothermal compressibility 
B(t) = transient compressibility 

c — phase velocity (frequency dependent) 
c.s = adiabatic (isentropic) phase velocity 
Cp = specific heat at constant pressure 
Cp— specific heat at constant volume 
G* = complex shear modulus 
Hv = retardation spectrum for volume deformation 

i = imaginary unit 
J n == Bessel Function 
K* =K' + iK" = complex bulk modulus 
/W* =M' + iM" = complex longitudinal wave modulus 
p = s + ia) = complex angular frequency 
Po = 5o + io)o = singular point of B* (p ) 
po = So — ia){) = complex conjugate of p () 
5 = Re p 
5 = Re po 
t = time 

u = displacement 
x = position coordinate 
a = Re r 
(3 = (i)lc = Im r 

r = a -\-ifi= complex propagation function 

y = CplCv = BrlBs = ratio of specific heats, or compressibilities 
tj = shear viscosity 
k = thermal conductivity 
p = density 

r = time, relaxation time 
v = dilatational viscosity 
co = Im p = angular frequency 

&>o = Imp 
93*(p) = normalized dispersion function pertaining to complex compressibility 
®(t)-«.£f-i[©*(p)/p] 

J£ = Laplace transform 
oS^ _1 = inverse Laplace transform 
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As a result of the complicated algebraic form of the apparent steady-state solutions and the diffi- 
culties which arise in obtaining the general solution [9] of the differential equation, the invalid 
assumptions mentioned in the last paragraph have been widely accepted. In many cases con- 
clusions were based largely on qualitative considerations (for example, references [10-13]) and 
approximations which facilitate simple calculations, but which shroud other important information. 
Recently, the analysis has been facilitated by the advent of fast digital computers, which make it 
a practical matter to trace the branches of the exact propagation functions over a wide frequency 
range for as many different cases as necessary to illustrate the complete behavior. 

2. Steady-State Functions 

2.1. Propagation Functions 

In this work we are principally concerned with the forced vibration (sinusoidal) solutions of 
the differential equation [4] (see appendix), 
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where u— u(x, t) is the displacement 2 of an element in the x direction at time t, c s is the adiabatic 
(or isentropic) phase velocity, k is the thermal conductivity, p is the density, taken here in the 
absence of any disturbance within the medium, Cy is the specific heat at constant volume, and y 
is the ratio of specific heats (Cr/Cy). c s -, K, p, CY, and y are taken to be constants of the fluid 
independent of x and t, which implies that these quantities are independent of the small fluctuations 
appropriate to the linearized equation (eq (1)). Equation (1) describes the infinitesimal (i.e., linear), 
uniaxial response of a hypothetical continuum of compressible, inviscid, heat-conducting fluid. 
This is probably the simplest case for which the role of thermal conductivity is considered in 
acoustic propagation. For our purposes eq (1) may be taken axiomatically and is used as a point 
of departure for subsequent arguments. 

We are concerned with obtaining the exact solution 3 of eq (1) subjected to the boundary 
condition w(0, t) = u () U(t) cos cot (along with other obvious statements with respect to initial 
conditions and behavior as x — > ») where U(t) is the unit step function. Although obtaining the 
Laplace transform, J?f[a(je, t)]= u(x, p), is straightforward procedure, obtaining its inverse 
u(x, t) appears to be extremely difficult. However, sufficient information can be obtained from 
u(x, p) to reveal that the exact solution is unbounded, or unstable, when 1 < y < 2. 

In view of the difficulties in obtaining the exact solution, we will assume the form of the 
steady-state solution of eq (1). From this we will obtain the frequency-dependent propagation func- 
tion I(p) and its equivalent response function Z?*(p), which appears in simpler form. The stability 
will be evaluated by applying Laplace transform theory to the response function, and verified later 
by obtaining the corresponding time-dependent (transient) function. 

For waves traveling in the positive x direction only with sinusoidal time dependence, the 
apparent 4 "steady-state" solutions of eq (1) may be expressed in the following form: 

u(x, t)=Re [Ae*> t - rx ] 8 =o = Ae- aa: cos (cat — fix) (2) 

where p = s-\-ico is the complex angular frequency, 5 and T=a-\-i/3 is the complex (frequency 
dependent) propagation function with /3 = cole, where c is the frequency dependent phase velocity. 



2 The equation for velocity at a point, v = v(x, t) . has identical form. 

:1 This type (including transient it-spouse) <>i analysis with respect to acoustic propagation is, in general, fairly recent in the literature. For the exact solution 
pertaining to a viscous, nonconducting fluid, see, for example, reference [141. 

4 We speak of eq (2) as an "apparent" steady-state solution inasmuch as it appears to be a steady state and satisfies the differential equation (eq (1)). However, 
it remains to be seen if eq (2) is the exact solution as t~*°°. 

5 The complex frequency, i> = s + iu>. is used for the purpose of applying the Laplace transform in the subsequent discussions. It is convenient to use the form 
exp (pt — Vx) in lieu of exp {Tx—pt) in order to facilitate the subsequent use of this transform without modification. 
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A is a constant (for each branch), which is taken to be real without loss in generality. 

Substitution of eq (2) into eq (1) gives the following frequency dependence of the complex 
propagation function T(p): 
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where a = pc 2 s C\ /k, or alternatively, a = C\ -jkBs, remembering c|= 1/pBs where Bs is the adiabatic 
compressibility. Equation (3) may be identified as the propagation function for "steady-state" 
solution of the Kirchhoff-Langevin equation [8, 9] with zero viscosity, or equivalently, with an 
infinite viscous Reynolds number. 

For waves with diminishing intensity in the direction of propagation, as required for passive 
fluids, a and ft must have the same sign. This requirement is met for all four branches given by 
eq (3). The two branches for which a and (3 are both negative are ignored in this analysis because 
they correspond to backward traveling waves and contribute no additional information to the 
properties of the fluid. The two branches implied by the i-power term in eq (3) are commonly 
called the acoustic and thermal branches. As customary, the acoustic branch is taken to be the 
one for which the absorption, a, vanishes in the low-frequency limit in contrast to an unbounded 
value in the thermal branch. This choice is consistent with the acoustic branch propagating adi- 
abatically in the low-frequency limit. 

Expressions for a and /3 as functions of real frequency are obtained tediously, producing 
expressions containing eight radicals each, which, indeed, exceed the practical limit of compre- 
hension and space available here. However, the values of these functions with frequency were 
computed for several examples and are illustrated in section 2.4 along with their corresponding 
complex response functions. 

2.2. Complex Response Functions 

For mathematical convenience the propagation function obtained in the last section is con- 
verted to the complex response function which describes the response to alternating pressure of 
any small volume element of constant mass within the acoustic field. The information obtained in 
these two representations is equivalent. 

The longitudinal wave propagation function and its corresponding complex response functions 
are related [15] by 

M* = K*+|c* = pp 2 /r 2 (4) 

where M* = M' + iM" is commonly called the longitudinal wave modulus, which is appropriate 
to describe linear response for uniaxial strain with sinusoidal time dependence. K* and C* are 
the complex bulk and shear moduli, respectively. Since we are concerned with inviscid, compress- 
ible fluids, which support no shear stress, G* may be taken to be identically zero from which follows 

B* = P/pp 2 , (5) 

where B* = B' — iB" '=1IK* is the complex compressibility. Thus, B* and F 2 have the same fre- 
quency dependence except for a factor of p 2 . The alternative expressions" equivalent to eq (5) 
in terms of real and imaginary components at real frequencies are 

B' = (/3 2 -a 2 )loj 2 p 

(6) 
B" = 2ap/a) 2 p. 



6 From inspection of eq (6) it is clear that there are only two distinct branches of B*(p) from the four branches of V (p) given by eq (3), since a and (3 always 
have the same sign. 
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As commonly used in acoustics (and in this discussion) the complex compressibility relates 
the dilatation to pressure in both amplitude and phase for sinusoidally time dependent excitation. 
B*(p) is not a property of the fluid alone, as are the isothermal and adiabatic compressibilities 
for a viscous fluid, but may also depend on the boundary conditions. Other response functions, 
such as the characteristic impedance, or admittance, of the medium, could be analyzed in this 
manner. Since these response functions are equivalent, the conditions for stability obtained through 
analysis of any of them remain the same. The complex compressibility appears to be the most 
convenient function for studying this particular case. 

The identity, 
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where c and d may take on all real values, is useful for establishing the inverse relations to eq (6) 
and certain other useful relations obtained subsequently. In this paper the notation z 1/2 is used 
for the double-valued function. Expressions of the form \z are interpreted to be single-valued 
functions with a nonnegative real part. The other branch will be denoted by — VT (with a non- 
positive real part). This convention so far does not distinguish the branches of z l/2 when the value 
of z is on the negative real axis. For a reason which will become apparent in the treatment of the 
special case y — 2 (sec. 4.2), the branch cut is taken along the negative real axis and is assigned 
the polar argument the value = — rr (instead of the more customary choice S=tt). Thus, for 
the first branch (— tt ^ < tt) we write z ,/2 = Vz, which has a nonnegative real part, and for the 
second branch (tt ^ 6 < 3tt), z ll2 = — Vz, with a nonpositive real part. This choice has the conse- 
quence that V— 1 = — i instead of + i as usual. The form ± Vz or+ Vz is used according to the 
convention that the upper sign applies to the acoustic branch, whenever applicable. 
The inverse relations to eqs (6), 
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are readily obtained from eqs (5) and (7). 

From eqs (3) and (5) the explicit frequency dependence of the complex compressibility is 
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which, after taking p = ico (5 = 0), becomes 
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The following real and imaginary components of B*(co) = B' — iW are obtained by applying eq (7): 
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Equations (10) and (11) are in keeping with our convention given in the last paragraph and con- 
sistent with our previously stated choice (sec. 2.1) that the acoustic branch be the one for which 
the absorption vanish in the limit co^ 0. This choice is consistent with B*{0) = B s for the acoustic 
branch. The two expressions of £*(o>) and B' ((d) are shown here for the purpose of distinguishing 
between the acoustic (upper sign) and thermal (lower sign) branches. The reason for the sign 
reversal at 7 = 2 is found in the fact that the roles of the two branches are interchanged in a way 
which should become clearer during further discussions. Values of y less than unity are not con- 
sidered because these would require the existence of a negative compressibility or specific heat. 
It will be shown in the next section that for 1 < y < 2, the apparent steady-state solutions are 
misleading because the fluids are unstable in that the corresponding "transient" functions are 
regenerative in time. Thus, in excluding the unstable solutions, eqs (10a) and (11a) are applicable 
only to the special case where y= 1. 

From initial inspection eq (9) appears to have a simple pole at the origin in addition to the 
branch-point singularities. A binomial expansion of the radical reveals that the simple pole in the 
acoustic branch disappears by virtue of a cancellation in terms of alp. The presence of the simple 
pole in the thermal branch is consistent with the unbounded absorption in the low-frequency limit. 

2.3. Evaluation of Stability from the Location of Singularities 

It is commonly assumed in wave propagation that transient solutions decaying with time, are 
conjugate to steady-state solutions decaying in the direction of propagation. In this section it will 
be shown that this assumption is not correct in general and is violated for the response functions 
under consideration when 1 < y < 2. 

Since the general solution of eq (1) is extremely difficult and apparently has not been obtained, 
we have chosen an alternate approach, namely, to study the transient 7 response function B(t), 
where B(t) is the ratio of time-dependent dilatation, to step function pressure. B(t) corresponds 
to the response function B*(p) in frequency space. (See eq (9).) The complex and transient func- 
tions are related by 

B*(p)=pJ?[B(t)] 



with the corresponding inverse 



B(t)=3?-i[B*(p)/p], (12) 



where J? and Jf' 1 are the Laplace and inverse Laplace transform operators. We say a solution of 
eq (1) is stable if the corresponding transient response function is bounded for all times t after 
removal of terms linear in t. This requires B*(p) to have singularities only on, or to the left of, the 
imaginary p axis. (See, for example, [16, 17].) Terms of B(t) linear in t (corresponding to simple 
poles of B*(p) at the origin) are permitted, although they are unbounded in the limit as t^oc. 
These terms are analogous to those for steady-flow shear within a Newtonian fluid. In general, for 
compliances (resulting from step-function stress) of the form ct n , the only acceptable values of 
n for stable response are O^ra^l. n<0 would give a singularity at t = implying an infinite 
instantaneous deformation from the application of stress, n > 1 characterizes a deformation which 
would continue to increase after removal of the stress (and without bound in the limit t— * °°), as 



7 The response to step-function excitation is commonly called the transient response function, in view of the fact that this function usually vanishes at large 
times. In this case the so-called "transient" response function B{t) often grows with time and is unbounded in the limit as / — * ». This behavior is contrary to that 
implied by the ordinary definition of transient, i.e., existing over only a short time. 
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can be seen by applying the Boltzmann superposition principle. < n < 1 corresponds to a de- 
formation which would recover completely on removal of stress but would increase without limit as 
long as the stress is applied. n = and n=\ represent purely elastic and purely viscous response, 
respectively. 

In addition to the simple pole at the origin, branch point singularities appear at the values 



(7 — 2) a 2iavy—l 
p = ± 

y y 

Thus, as required for real time-dependent response functions, singularities off the real axis are 
always complex conjugates. The locus of these singularities with respect to 7 is a circle obeying 
the relation 

s 2 + co 2 o = a 2 . (13) 

This plot is shown on figure 1 in which the values of y are given for the intercepts. Singularities 




FIGURE 1. Locus of the branch-point singularities of the 
complex response function, B*(p), given by eq (9). 

Open circles corresponding to y — I and y^x indicate exclusions from the locus. 

are restricted to be on and to the left of the imaginary axis when y ^ 2. When y is exactly one, the 
radical in eq (9) is a complete square and the point at p = a is regular. However, as soon as y differs 
from unity (even infinitesimally), a pair of branch points appears and the corresponding instability 
sets in. Accordingly, the point y=l is an isolated stable point. At 7 = 00 the branch-point singu- 
larities become either simple poles or vanish entirely. (See sec. 4.3.) Accordingly, the points p = a 
(corresponding to 7=1) and p = — a (corresponding to 7 = 00) are excluded from the locus as indi- 
cated by the open circles in figure 1. Since the singularities are paired symmetrically about the 
real p axis, the transient response will be oscillatory. (See sec. 3.2.) When the singularities are 
paired about the negative real axis as with 7^2, the oscillations damp out with time. When the 
singularities are paired about the positive real axis, as they appear for 1 < 7 < 2, the oscillations 
amplify with time and are unbounded as t — > 0°. 

2.4. Graphical Display of the Functions 

The values of the apparent steady-state functions plotted in this section were obtained in 
three different ways for the purpose of checking the validity of the results. The functions were 
programmed in OMNITAB and also in FORTRAN using both the real and complex modes. With 
OMNITAB and FORTRAN using the real mode, eqs (11) were used to compute the values of the 
real and imaginary components of the normalized complex response functions. These were con- 
verted to the values of the normalized propagation functions using eqs (8) with j3 = co/c. With 
FORTRAN using the complex mode eqs (10) were used to compute the values of the normalized 
complex response function, which were converted to the values of the normalized complex propaga- 
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tion function using eq (5) rewritten in the form 



r(co)=-coVp£*(co). 

Since we have taken the branch cut of the function z 1/2 in a manner consistent with V— 1 =— i 
(see sec. 2.2.), whereas the FORTRAN complex mode assigns the value i to the operation V— 1, it 
was necessary to include a statement to make the corresponding inversion at y = 2 over the range 
cd > a (see sec. 4.2.). Figures 2 and 3 give the computed dimensionless complex response functions 
plotted against dimensionless frequency over the interesting range including the value corre- 
sponding to o) = a for different values of y. The principal distinction between these two figures 
is that figures 2 and 3 pertain to unstable and stable fluids, respectively. Figures 4 and 5 give the 
corresponding dimensionless propagation functions, which are, of course, more directly realizable 
experimentally and discussed in the literature. The functions for which the fluids are unstable 
are included because they are illustrations of functions which have been considered to be stable. 
Figures 2 and 3 with y= 1.5 are an example in the range of 1 < y < 2 over which the fluids are un- 
stable. The unstable fluids are the only ones (except when y = 2) which would seem to propagate 
monotonically and continuously between adiabatic and isothermal limits. Figures 3 and 5 with 
y=l, 2, and 2.5 are examples for which the fluids are stable. From these figures the nature of the 
coalescence and inversion of the branches at y = 2 (see sec. 4.2) is revealed more clearly. 

With the exception of the special case of y=l, the compressibility always vanishes in the 
high-frequency limit of the acoustic branch for stable fluids. This corresponds to so-called "anom- 
alous" propagation in that the phase velocity is unbounded in the high-frequency limit. Although 
the acoustic branch propagates adiabatically in the low-frequency limit, and the thermal branch, 
isothermally in the high-frequency limit, the functions are never continuous between these limits 
for stable fluids except when y = 2. In this case, the infinite derivatives at the value corresponding 
to a) = a in all of the functions are consistent with the existence of a branch-point singularity at 
this value. This singularity also accounts for the coalescence of B'(oj) and B"(qj) in the range 
co ^ a and oj ^ a, respectively, with y = 2. 

The dotted lines shown with the case for y = 2.5 were calculated for the components of the 
viscoelastic response function 

B*(iw)=-^— (14) 

1 + IOJT 

where r is the relaxation 8 time for a single element relaxation spectrum. The limiting values of 
eq (14) for small and large p are the same as those for eqs (lib) and (lie) using the upper signs. A 
single relaxation manifests the narrowest dispersion region possible with respect to any non- 
negative relaxation distribution. Accordingly, it is clear from inspection alone of figures 3c and 5c 
that the in viscid functions under discussion cannot be represented by any distribution of relaxation 
(or retardation) times in the form 

™-W."f^* (i5) 

with Hv(t) nonnegative. For a single relaxation time Hv(r) is a delta function multiplied by a con- 
stant. Since eq (15) applies only to functions with singularities on the negative real p axis, it is not 
applicable, in general, to the functions under discussion. 



8 This usage is apparently common in acoustic nomenclature. According to viscoelastic nomenclature [18], this quantity should be called a retardation time when 
used in connection with compliances as with eqs (14) and (15) in contrast to a relaxation time, with moduli. 
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FIGURE 2. Real and imaginary parts of the normalized 
complex response functions as functions of normalized 
frequency with y = 1.5 as an example in the range of 
1 < y < 2 for which the propagation is unstable. 
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4. Normalized propagation functions corresponding 
the functions in figure 2 (unstable propagation). 
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FIGURE 3. Real and imaginary parts of the normalized 
complex response functions as functions of normalized fre- 
quency with (a) y=l, (b) y = 2, and (c) y = 2.5 as examples 
at y = 1 , and in the range ofy 2= 2 for which the propagation 
is stable. 
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FIGURE 5. Normalized propagation functions corresponding 
to the functions in figure 3 (stable propagation). 



3. Transient Response Functions 9 

3.1. General Integral Formulation 

In the previous paragraphs we have given a necessary and sufficient criterion to distinguish 
between stable and unstable fluids. For the purpose of obtaining additional insight into the detailed 
behavior of the solutions, it is useful to study more closely the transient response functions which 
can be written down explicitly as inverse Laplace transforms of B*(p)lp. It is convenient to split 
the expression for B*(p)/p (see eq (9)) as follows: 



B*(p) _ yB s 

P 2 



l + a + mEl 



where 



33* (p) = 1 
P P 



yp p 



(16) 

(17) 



Application of the inverse Laplace transform to eq (16) yields: 



^{1)=^ [l + ar+ 23(f)]. 



(18) 



This reduces the problem to finding the explicit time dependence of 93(0? a function which is 
double valued, owing to the V2-power in 53* (p). 



If we put 
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we may write eq (17) as follows: 



po = s + ico {) 
po = So — icoo, 



^^ = -A(p~ Po) (P - Po) Y* = -3 [p 2 - 2s»p + a 2 ] «\ 



P P 

which can be decomposed still further into 



(19) 



«M. flW .2s B(p)+ ^ (p) , 



where 



D{p) = [p 2 -25op + a 2 ]-'/ 2 . 
From equation 5.3.34 in [19] we find 

^-'{£>(p)}=e«,/o(a>„ C ). 



9 See footnote 7. 
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Writing ± to allow for the double-valuedness of the V2-power function, we obtain 

±93(0 =e*7o(<»o*) - 2*o J e s » T Jo(a) {) T)dr 



+ 



^Wo^ 



(ojot) drdt' . 



Partial integration in the double integral and application of equation 11.3.5 in [20] yields the final 
result: 

± 93(t) = e *[(! -Sot)Ma> t) ~ (OotJduot) ] 



J o 



+ (a 2 f-s«) e 5 " T Jo(o) {) r)dT. 



(20) 



Remembering our convention that the upper sign pertains to the acoustic branch, the exact 
form of the transient response function B(t) for y ^ 2(s<> ^ 0) is obtained from substituting eq (20) 
into eq (18): 

B(t) =^2^|l +aiT e*"'(l - s t)Ma>ot) - Q>otJi(a>ot) 

+ (a 2 *-s ) J e«*Ma>oT)dT ], y & 2. (21) 

When 1 ^ y < 2 (s„ ^ 0), eq (21) is valid also except T becomes ±. 



3.2. Approximate Behavior for Large t 

Although eq (21) is the general relation applicable for all values of t and y ^ 1, it is difficult to 
obtain comprehensible information with respect to the explicit asymptotic behavior as t —* o° 
which is of particular interest here. To obtain this information we take the expansions of B*(p)lp 
about the three singularities p = 0, p = po, and p = po. (From the general theory of Laplace trans- 
forms [21] the asymptotic behavior of y(t) for large t may be readily determined by the behavior of 
the leading terms of J?[y(t)] =y*(p)lp expanded about the singularities in the p plane.) 

In this connection we are interested in evaluating eq (16) (sec. 3.1) which is 



where 



B*(p) yfi.s 
P 2 


ri + a + 93*(p)" 

-P p 2 p . 


33*(P)_1[ 1 

P PL 


, 2(y-2) , an 

yp p 1 . 



1/2 



(16) 



(17) 



Since the inverse transforms of the first two terms in the brackets of eq (16) are readily available, 
we are only concerned with the appropriate expansions of eq (17). 

a. Expansion about p = 

For the expansion about p = 0, we rewrite eq (17) as follows: 

s #*(p) _ a ( l , 2(y-2)p | g^"» 



ay 
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For small p the above equation may be expanded by the binomial theorem indicating double 
valuedness by the + sign. 



93*(p) __^a , y-2 



P 



p A yp 



■> 



The dots indicate regular terms which do not contribute to the asymptotic expansions of either 
B(t) or 93 (t). The corresponding expression for 93 (^) , for large t, is then 

where the subscript indicates that only the contribution from p = is included. 



b. Expansion about p = po 

In this case we rewrite eq (17) as 

93*(p) a r/ w xni , 

—y-^ = ~2 [(p-po)(p-po)] 1/2 

which when rewritten in appropriate form for expansion about p = po becomes 



»*(p) 
P 



a(po — po) 1 ' 2 \ go-go 



1-f 



P~ Po 



1/2 



P5 



1 + 



P-Po V 
Po / 



(P"Po) 1/2 , 



where all the terms within the brackets are regular at p = po. Consequently, 

SB*(p)_a(po-po)' /2 



PS 



(p-po) 1/2 (l + . • •), 



where the omitted terms in the right-hand parenthesis are ignored because they give rise to a 
series expansion of 93(f) in decreasing powers of t, thus contributing less as £— >°°. From reference 
[21] we find 



8i(t) 



a jpo — 

2pl V 77 



Po 



e p t t -3l2^ t^OO 



as the leading term for the contribution from p = p . 

c. Expansion about p = po 

The expansion about p =po is obtained from the above expression for 93i(0 by interchanging 
Po and po, which is tantamount to taking the complex conjugate of 93i(f). Thus, 



2Po > 77 



e~ pot t~ sf2 , tr-*a° 



pertains to p = po. 
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d. Summation of Contributions 



Adding the above contributions (93o, 85,, and 93,) together along with the contributions 
from the two left-hand terms in eq (16), we obtain the following general (y ^ 1) asymptotic expres- 
sion for B(t) as t —* <*>: 



B{t) 



yB. 



'-{ 



1 + at 



at +■ 



t -3,2 Re iXE»EE eP « 



Vtt 



Po 



Alternatively, the following equivalent expression may be obtained in terms of real quantities: 



B(t) -^1 
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at + z h 



(55 + W5) 2 V 7T 
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(22) 



where the upper sign now pertains to the acoustic branch for all y ^ 1. 
For the stable class y ^ 2 the asymptotic behavior is 



*M 



yB s 



1 + at 



at 



7-2 
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Asymptotes linear in / (although unbounded in the limit t— * °°) are consistent with stable solutions 
for the reason given in section 2.3. At large t sinusoidal oscillations about these asymptotes are 
damped by the factor t~* ,2 e s>,t for which s ^ 0. For the unstable class 1 < y < 2 these oscillations 
regenerate by the same factor for which s > 0. 

4. Special Cases for Different Values of y 

The following discussion includes the special cases of the functions for different values or 
ranges of y. In these cases the functions may reduce to considerably simpler form and may also 
reveal distinct characteristics not immediately obvious from the more general form. 

4.1. Stable Case y = l 

When y= 1, the reduced expression for eq (9) becomes 



B*(p) 



B s 



1-f- 



1-2 



,'y=i. 



(23) 



Thus, it is apparent that B' = Bs and B" = for the acoustic branch, and B'=0 and B" = C\/kco 
for the thermal branch. These values are consistent with a constant adiabatic (or isothermal) 
phase velocity with no absorption in the acoustic branch. The thermal branch propagates in a 
manner analogous to shear propagation within a Newtonian fluid [2] in that 



a = (3-- 



copC\ 



2k 



The steady-state functions for this case are displayed in figures 3a and 5a. 

At y=l the branch point singularities vanish because the square in eq (19) is completed. 
Accordingly, this point is an isolated stable point and the response is characterized by the simple 
pole at the origin. 
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The following expression for B(t) (which is exact for all t) is obtained by operating on eq (23) 
directly: 



B(t) 



Bs 
2 



[I -\- at ± (I -at)], y=l. 



It is apparent from the above equation that when y=l, the thermal branch (lower sign) of B(t) 
is analogous to the shear compliance of a Newtonian fluid subjected to a constant stress. Transient 
terms linear in t (although unbounded in the limit £— »o°) are consistent with stable solutions for 
the reason given in section 2.3. 

4.2. Stable Case y = 2 

When y = 2, eqs (lib) and (lie) reduce to 



B'(oj)=Bs\ 
B"(cj)=B, 



VF? 



-h 



1 



-fcM- 



,=2. 



It is apparent from the above equations that B' is single valued and constant over the range cj ^ a, 
and B" is single valued over the range co ^ a, in a manner consistent with the branch point singu- 
larities at p = ± ia. In this case the two branches of B' coalesce for ^ co ^ a by virtue of the fact 
that (y— l)B s = B s . The steady-state functions for this case are displayed in figures 3b and 5b. 
Since this case corresponds to d=0 in eq (7), and because of the single-valuedness mentioned 
above we have some liberty in designating the branches with respect to acoustic or thermal over 
the range co > a. If we had used the customary convention in setting V— 1= z, in lieu of V— 1 =— i 
(see discussion after eq (7), sec. 2.2), the branches (as shown on fig. 3b) would be interchanged 
over the range co > a. Since this case has been shown to be stable, we have chosen to deviate from 
the customary convention. This results in the limiting behavior of B*(^) being the same as that 
for the stable class y > 2 in that 



Acoustic 



£*(oc)=0, y^2 



1 <y<2. 



Thermal B*(™)=yB s , 

The response function B{t) valid for all t (expressed in integral form), and approximate for 
large t, is obtained by setting so = in eqs (21) and (22), respectively. The asymptotes are 

lim B(t)=B s [l + at + at], y = 0. 

At large t the sinusoidal oscillations about these asymptotes are damped by the factor t~' sl2 . This 
behavior, which is characterized by the branch points, is distinct from that of undamped resonance, 
characterized by simple poles also paired symmetrically on the imaginary p axis. 

4.3. Stable Case y = oc 

Since there is no phenomenological bound on the maximum value of y, we are considering 
the role of all possible values for which 1 ^ y ^ oo. (For example, y is apparently singular at the 
critical point [22] in that Cp is more singular than Cy or alternatively, B r is more singular than Bs.) 
The general form of B*(p) taken from eq (9) in the limit as y — > °° is 



lim B*{ P )=^ 



1+- 



1 + 



01 



>+a)/]} 



y(p+a) 



(24) 
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where the O symbol is the one used in asymptotic theory. The value of B*(/>) at y = oo depends on 
the dependence of B s on y in the limit as y — > ». The behavior of B*(p) at y= oc is summarized in 
table 1, where the acoustic and thermal branches are treated separately. It is interesting to observe 
from comparison of eqs (9) and (24) that the branch-point singularities coalesce into a simple pole 
at p = — a corresponding to y = °°, and, as indicated by table 1, the simple pole at p = — a vanishes 
for finite values of Br (fis = 0). 

The general form of B(t) in the limit as y— > oo corresponding to eq (24) is 



B T 



lim B(t)=- n 



ll + atT l + at-o(-(l-e- at )) j- 



The behavior of B(t) at y = °°, is also summarized in table 1. Although B(t) may be unbounded at 
y = oc for all values of t under certain conditions as shown, this case is considered to be a stable 
one because the response does not grow with time after removal of the terms linear in time. A 
medium of infinite compressibility, for example, a vacuum, does not, necessarily, have intrinsic 
energy sources as implied by our criterion for an unstable fluid. 

Table 1. Behavior of B*(p) and B(t) at y = °° 
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y~ *x 


B*(p), Acoustic 


B*(p), Thermal 
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fl. s = (Constr\ 
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B s (\ -e~'") 
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BtO + at) 

00 

oc 



5. Summary of Asymptotic Values and Relations Between Asymptotics of B*(p) and B(t) 

Table 2 summarizes the asymptotics for B*(p) and B(t). The limiting values for B*{p) with 
/; = 0, and p— > o° are taken from the appropriate binomial expansions of eq (9). The values of B(t) 
with £=0, and t— > oo are taken from eqs (21) and (22), respectively. The ±(°°) notation indicates 
regenerative oscillations modulated symmetrically about the value zero with the amplitudes going 
to infinity as t —> oc. 

The following give the relations between the asymptotes of B*(p) and B(t) 



B*(oo) = B(0), 
fi*(0) = fl(oo), 
B*(0) # 5(oo). 



y ^ 2, y = 1 
1 < y < 2 



The first is readily obtained from Laplace transform theory without recourse to any relations 
given in this paper. The first two are always obeyed in passive response functions. The reason the 
second is not obeyed in the unstable class (1 < y < 2) is simply because the rightmost singularities 
appear to the right of the imaginary axis of the p plane. 





Table 2. Summary of asymptotics o/B*(p) and B(t) 




Stable class 
(y=l, y^2) 


Unstable class 
(1 <y<2) 




Acoustic 


Thermal 


Acoustic 


Thermal 


B*(0) 
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B(0) 
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Bs(y-l+ya/p) 

B s (y- 1 +yat) 
B T 
B, 


B s 

B s ±(oo) 

B T 
B t 


Bs(y-l-rya/p) 

B s (y-\ +yat)±(?c) 





For meaning of ± (°°), see text. 
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6. Concluding Remarks 

The analysis of the differential equation for acoustic propagation in a continuum of an inviscid, 
compressible, heat-conducting fluid yields solutions which display unexpected properties. Although 
the "apparent" steady-state, forced-vibration solutions correspond to a physically reasonable 
behavior in that they predict oscillations which always decay in the direction of propagation, the 
corresponding "transient" response functions predict a response (to step-function excitation) 
which at long times may oscillate with increasing amplitudes. A regenerative "transient" response 
implies that the exact solution of the wave equation oscillates with increasing amplitude which is 
unbounded in the limit of infinite time, and, accordingly, a true steady-state solution does not exist. 
Since we restrict "acceptable" solutions to those which predict physically realistic behavior for 
passive, or stable, fluids, we reject the regenerative, or unstable^ solutions. 

The following statements are considered by the authors to be the most important consequences 
of the foregoing analysis: 

1. Contrary to the popular concept, the phase velocity can never vary continuously with 
frequency from adiabatic to isothermal values in a compressible, inviscid, heat-conducting fluid, 
except when the ratio of specific heats is exactly equal to two. 

2. Propagation within the above fluid is stable only when either of the following conditions hold: 

(i) y=l 

or 

(ii) 7^2. 

With the unstable class (1 < y < 2) the "steady-state" phase velocity of the acoustic branch 
appears to vary continuously with frequency from adiabatic to isothermal values. With this class 
the apparent steady-state solution is actually conjugate to a regenerative "transient" solution 
unbounded in the limit as t — > ». 

With fluids where the ratio of specific heats lies between the values one and two, stability may 
be achieved by the inclusion of a minimum viscosity. In a previous paper [1] the stability criteria 
were developed in terms of thermophysical properties applicable to all values of y in this more 
general case. These requirements are summarized as follows: 

7=1 
or 

(2— y)* 



v+ h) 



1, y>\ 



where v and r\ are the bulk and shear viscosities. Either of the above conditions is sufficient to 
insure stability, and these two are the only stable cases. 

Although the inviscid fluid is a hypothetical one (since all fluids presumably exhibit a finite 
viscosity), the unstable fluids cannot be stabilized by the inclusion of a viscosity of infinitesimal 
magnitude. The results from reference [1] clearly show that a lower bound exists on viscosity to 
insure stability when 1 < y < 2. The above stability requirements were found to be obeyed for all 
the fluids checked except the molten metals, including mercury, and possibly the superfluids. 
These results indicate that in some cases different assumptions would have to be made in deriving 
the propagation equation (Kirchhoff-Langevin equation) for compressible, viscous, heat-conducting 
fluids. 

In general an alternative evaluation of stability may be done by an appropriate qualitative anal- 
ysis using the Kronig- Kramers [23], or the equivalent Ginzberg [24], relations, which relate the real and 
imaginary parts of passive complex functions. (The Ginzberg relations apply specifically to velocity 
and absorption.) In this manner it is often possible to relate the sign of one complex component 
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to the sign of the trend of the other. For example, referring to figures 2 and 3c, it is possible to show 
that B' (acoustic) must be monotonically decreasing, and B' (thermal), monotonically increasing, 
with frequency for stable propagation. These trends are not obeyed in figure 2 because correspond- 
ing propagation would be unstable. The results with respect to stability are equivalent to those from 
the foregoing analysis because the Kronig-Kramers relations imply the existence of all singularities 
on, or to the left of, the imaginary axis of the complex p plane. 

The results of this work and those from the more general case [1] indicate the possibility that 
the apparent steady-state relations obtained from either phenomenological or molecular theories 
may be inconsistent with stable propagation. In particular, whenever thermoconductivity is included 
in a model, stability may be questionable. As illustrated here, the fact that the apparent steady- 
state solutions predict diminishing intensity in the direction of propagation is necessary, but not 
sufficient, for stable propagation. 

7. Appendix. Derivation of the Propagation Equation 

Equation (1) of the text may be obtained from the hydrodynamic equations [25] applicable to a 
Newtonian fluid under very general conditions. These equations are the equation of continuity, 
the Navier-Stokes equation, and the Fourier-Kirchhoff-Neumann energy equation. After dropping 
viscosity terms and linearizing (replacing DjDt by d/dt, dropping products of derivatives, and 
taking p constant when it appears in a product with a quantity which goes to zero with vainishing 
motion), these equations for motion in the x direction are 

dp dv rt 

— + p— =0 (A 1) 

dt H dx v } 

dv dP 

'*—■£ < A - 2 > 

p ^ P f+™ (A . 3) 

dt dx dx 1 

where v, 6, P, and e are the velocity, temperature, pressure, and internal energy per unit mass, 
respectively. 

For an inviscid fluid we take e = e(p, P) and 0=0(p, P) to be single-valued functions of the 
indicated arguments. From eqs (A. 2) and (A. 3) the following ordinary thermodynamic relations are 
implied: 



(■ 



de\ 
dPJ, = 


yBsCy 
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p 1 pap 



Expanding eq (A.3) by the chain rule, we obtain 

(de\ dp /de\ dP = _ p dv, (d0\ d^p (dd\ d^P 
P \dp)pdt P \dP) p dt dx K \dp)pdx* K \dP) p dx> 

in which the derivatives operating on e are removed by substituting the above thermodynamic 
relations, and those on 0, by substituting the definitions 



Bs ~ l P w)s~^w)e' a "-~pO; 



yp 
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Differentiation of this result with respect to x and t and substitution of eqs (A.l) and (A.2) to replace 
dP/dx and dp/dt gives 



dt 



d 2 v 

dt 2 ' 



d 2 v 

dx* 



pC\ dx 2 



d 2 v c 2 s d 2 vl 
dt 2 y dx 2 \ 



(A.4) 



Equation (1) of the text may be obtained by substituting du/dt for v. One time derivative may 
be integrated out by taking a definite integral from — oo to t and setting u(x, —°°) and all of its 
derivatives equal to zero. The result (eq (1) of text) is simply a substitution of u for v in eq (A.4). 
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